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Abstract 

This paper addresses robust communication on a fading relay channel in which the relay is connected 
to the decoder via an out-of-band digital link of limited capacity. Both the source-to-relay and the source- 
to-destination links are subject to fading gains, which are generally unknown to the encoder prior to 
transmission. To overcome this impairment, a hybrid automatic retransmission request (HARQ) protocol 
is combined with multi-layer broadcast transmission, thus allowing for variable-rate decoding. Moreover, 
motivated by cloud radio access network applications, the relay operation is limited to compress-and- 
forward. The aim is maximizing the throughput performance as measured by the average number of 
successfully received bits per channel use, under either long-term static channel (LTSC) or short-term 
static channel (STSC) models. In order to opportunistically leverage better channel states based on 
the HARQ feedback from the decoder, an adaptive compression strategy at the relay is also proposed. 
Numerical results confirm the effectiveness of the proposed strategies. 

I. Introduction 

Consider the fading relay channel model shown in Fig. \T\ in which an encoder communicates to 
a decoder through a relay that is connected to the decoder via an out-of-band capacity-constrained 
backhaul link. Both the source-to-relay and the source-to-destination links are subject to fading. 
The motivation for this model comes from the uplink of cloud radio access networks [[DED, in 
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which base stations (BSs) operate as soft relays that communicate with a central decoder via 
a digital backhaul links. In this scenario, the central decoder performs decoding based on the 
compressed signals collected from the connected BSs. With regards to this application, in the 
model of Fig. [TJ the encoder represents a mobile station (MS); the relay is the BS in the same 
cell, which is connected to the (central) decoder via backhaul link; and the signal Y represents 
the compressed signals collected by the decoder from BSs belonging to other cells. The signal 
Y can be seen as side information available at the decoder when designing the encoder (i.e., the 
MS) and relay (i.e., the BS). 

The fading relay channel was investigated in and flU in the absence of the direct link 
between the source and the destination and assuming fading also on the relay-to-destination link. 
In J3]|, various relaying protocols including decode-and-forward, quantize-and-forward and hybrid 
amplify-quantize-and-forward were developed in combination with layered broadcast coding 
(BC). This work was extended in 01 by studying infinite-layer coding at both the source and 
the relay in conjunction with decode-and-forward relaying. 

The fading relay channels with a direct link between the source and the destination (as in 
Fig. CO) was studied in [[5||-[[T2|. The works in [El-CD solved the problem of optimizing the 
compression strategy at the relay under the assumption of perfect channel state information 
for multi-antenna terminals. In the presence of uncertainty on the fading coefficient S, layered 
approaches that adopt a competitive, rather than average, optimality criterion are derived in 
|fTT| and 021 assuming no hybrid automatic retransmission request (HARQ). In all the previous 
works, the feedback link in Fig. \T\ was not included. This link is used in this paper to enable 
HARQ. 

A. Contributions 

In this work, motivated by cloud radio access applications as mentioned above, we study 
the system in Fig. \T\ assuming that the relay performs compress-and-forward. We propose to 
combine two key strategies to mitigate the impact of the fading on the source-to-relay and source- 
to-destination links, namely, HARQ and BC. With HARQ, the decoder requests retransmission by 
sending feedback information to the encoder and the relay regarding the outcome of the decoding 
process. Specifically, the incremental redundancy HARQ (IR-HARQ) consists of the transmission 
of additional parity bits in case of failed decoding [13]. With BC [[T4 l |- lfT6l . instead, one allows 
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for variable-rate decoding that opportunistically adapts to the actual fading state conditions. 

Multi-layer HARQ strategies have thus the advantage of allowing for variable-length transmis- 
sion and variable-rate decoding, and were introduced in [fPTl for point-to-point fading channels. 
As in IfTTll . we aim at maximizing the average throughput and distinguish two scenarios, namely 
short-term static channel (STSC) and long-term static channel (LTSC). Moreover, for the LTSC 
scenario, we propose an adaptive compression method at the relay that is able to opportunistically 
leverage better fading state based on the feedback information received from the decoder. The 
effectiveness of the proposed multi-layer HARQ strategies is confirmed via extensive numerical 
results. 

The paper is organized as follows. We state the system model in Sec. [IT] and establish the 
problem formulation in Sec. Hill After describing the proposed multi-layer HARQ strategies with 
a constant compression gain and adaptive compression gain in Sec. [IV] and Sec. [Vl respectively, 
for the LTSC model, we extend the discussion to the STSC model in Sec. ED Numerical results 
are provided in Sec. IVIII to demonstrate the performance gain of the proposed multi-layer HARQ 
strategies. 

Notation: We adopt standard information-theoretic definitions for the mutual information 
I(X; Y) between the random variables X and Y, and conditional mutual information I(X; Y\Z) 
between X and Y conditioned on random variable Z ifLSl . All logarithms are in base two unless 
specified. We use Ex[-] to denote the expectation over X. For a real number x, we define a 
function [x] + = maxjx, 0}. 



We consider the fading relay channel depicted in Fig. [TJ in which the relay is connected to the 
decoder via a digital link of capacity C max . In order to enable HARQ, after each transmission 
block (or slot), the decoder sends feedback information to the encoder acknowledging, or not, 
successful decoding. This feedback link is assumed to be error-free. 

A. Channel Model 

The signal V t)i received by the relay in the ith symbol, i — 1, . . . , n, of the tth transmission 
slot is given as 



II. System Model 




(1) 
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Figure 1. A fading relay channel with out-of-band relay-to-decoder communication and a feedback link to enable HARQ. 



for t = 1, . . . , T, where D t is the fading coefficient in the tth time slot, X t>i ~ A/"(0, P) represents 
the signal transmitted by the encoder, E t>i ~ jV(0, 1) is the additive noise at the relay, and T 
is the maximum tolerable delay for the HARQ process. As it will be detailed, upon correct 
decoding at the destination, the HARQ process is stopped, and T is the maximum number of 
overall transmissions allowed for the same data packet. We assume that the block size n is large 
enough to enable the use of information-theoretic limits. The notation D t has been chosen with 
reference to the cloud radio access application in which D t represents the direct channel to the 
local BS. 

The symbol Y t)i received by the decoder in the ith symbol, % = 1, . . . , n, of the tth transmission 
slot is 

Yt,i — y~StX ti + Z t ,i, (2) 

for t — 1, . . . , T, where S t is the fading coefficient in the tth time slot, and Z ti ~ A/"(0, 1) is 
the additive noise. The notation S t is a reminder that in the cloud radio access application, S t 
represents the side information channel (see Sec. H). From now on, we omit the symbol index i 
for notational brevity. 

Following IfTTTl . depending on the channel coherence time, we distinguish two scenarios: i) 
short-term static channel (STSC); and ii) long-term static channel (LTSC). With LTSC, the 
channels remain fixed over all the, at most, T transmission blocks used for the current data 
packet, that is, 

D t = D and S t = S, for aU t = 1, . . . , T. (3) 
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In contrast, with STSC, the channel changes independently from block to block. We first study 
the LTSC model ® in Sec. |IV] and Sec. El and then consider the STSC case in Sec. EH We 
assume that the fading coefficients D t in (OQ) and St in © are independent, and have arbitrary 
CDFs F D (d) = Pr[A < d] and F s (s) = Pr[S t < s] with finite powers p D = E[\D\ 2 ] and 
Ps = E[|5| 2 ], respectively. 

The realization of the fading coefficients St is known only to the decoder, while that of the 
fading coefficients D t is available at the decoder as well as the relay. In order to study the 
effect of the local CSI at the encoder, we will consider both cases where the encoder knows the 
realization of the "direct" fading channel D t to the relay, e.g., through feedback, or not. 



B. Relay Operation 

The relay compresses its received signal V t and sends a description W t to the decoder. Without 
claim of optimality, we assume a Gaussian test channel (see, e.g., [18]) as 



W t = Jd~ t V t + N t , (4) 



where a t is a non-negative compression gain and N t ~ A/"(0, 1) represents the compression noise. 
Using binning for distributed source coding at the relay by leveraging the side information © 
at the decoder, the latter can recover the description W t as long as the inequality 

I(V t] W t \Y t ) = ^\og(l + a t (^ + T ^^)^ <C max (5) 

is satisfied [18, Ch. 11]. Due to the mentioned CSI limitation, the relay should compute the 

compression gain a t as a function of the realization of the local fading D t without having 

information about the fading state S t in ©. Therefore, in order to guarantee that the decoder 

can always recover W t regardless of the realization of the fading coefficient S t on the side 

information ©, one needs to set the compression gain a t so that ([5]) is satisfied even for the 

minimum value s mm in the support of Fs(s) (i.e., s mm = inf{s : Fs(s) > 0}). This leads to 

(3(1 + s min P) a _ 
l/d + (1 + s mm /d)P 

where = 2 2Cmax — 1, for D t — d. We will consider different strategies for the choice of the 
compression gain a t in Sec. [TV] and Sec. El 
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C. Multi-Layer Hybrid-ARQ 

Following IfPTl , the encoder uses a two-layer BC transmission strategy coupled with HARQ, 
which is described next. The encoder wishes to deliver two messages Mi G {1, . . . , 2 nRl } and 
M 2 G {1, . . . , 2 nR2 }, which are independent and uniformly distributed, to the decoder. To this end, 
it maps message M; to a n-symbol codeword X™ t for 1 — 1,2. We assume independent Gaussian 
codebooks across the T blocks, that is, the codewords X™ t are independently generated with 
i.i.d. symbols A/"(0, 1) for all blocks t = 1, . . . , T. To describe the multi-layer HARQ strategy, 
we distinguish the following two transmission modes: i) BC mode; and ii) single-layer (SL) 
mode. In the BC mode, the encoder transmits the superposition 



'aPX u + VaPX 2 >, 



(7) 



for each symbol, where a G [0, 1] and a = 1 — a represent the fractions of powers allocated to 
the first and second layers, respectively. In contrast, in the SL mode, the encoder transmits only 
the second-layer codeword with full power P, and the transmitted signal X t is written as 



X t 



PXo t . 



(8) 



In the first slot t = 1, the encoder emits the signal X% in the BC mode © and the relay sends 
the compressed version W\ in © of the received signal V\ to the decoder. At the completion of 
the slot, the decoder first tries to decode the message Mf, if successful, it cancels the codeword 
Xi i from the received signal and attempts to decode message M 2 . Decoding is based on the 
received signal in slot t, which can be written as 

V a tDt 



Y, 



W t 
Y t 



X t + 



Z t 



(9) 



for t = 1,...,T. 

The decoder informs the encoder and the relay about the number of layers that were correctly 
decoded. If both messages are not correctly decoded, in the next slot t = 2, the encoder sends 
incremental redundancy information for both layers using the BC mode 0. Note that incremental 
redundancy entails that, as mentioned, the codebooks used at different blocks are independent 
(see, e.g., lfT3l0 . Instead, if only the first layer M 1 was decoded in the first slot, the encoder 
transmits incremental redundancy information only for the second layer by using the SL mode 
d8]). This process lasts until either both messages Mi and M 2 are decoded successfully or the 
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maximum number T of transmissions is reached. Therefore, if a message Mi is not decoded 
until the Tth slot, outage is declared for layer I. 

III. Problem Definition 

The problem of interest is the maximization of the expected throughput 77 as measured by the 
average number of successfully received bits per channel use. Using the renewal theorem (see, 
e-g-> 03 ), we can calculate the expected throughput 77 as 

E |R] 

V = ^rr, (10) 
E [L] 

where E [Ft] is the average rate decoded in a HARQ session, which consists of at most T 
transmissions, and E [L] is the expected number of transmission blocks for HARQ session. 
Expectations are taken with respect to the fading coefficients D t and S t . These quantities can 
be computed as lfT7ll 

E [R] = R x (1 - pl ut (T)) +R 2 (l- p 2 ont (T)) , (11) 

T-l 

and E [L] = £ tp 2 dec (t) + T {p 2 dec (T) + p 2 ont (T)) , (12) 
t=\ 

where the probabilities p l out (k) and p l dec (k) are defined as 

p l ont (k) = Pr [Mi is not decoded until slot k] , (13) 
and p l dcc (k) = Pr [M t is decoded in slot k) . (14) 

The probabilities p l out (k) and p l dec (k) depend on the parameters R\, R 2 and a as will be clarified 
in the next sections. The problem of maximizing the average throughput 77 is then formulated as 

maximize r](R 1 , R 2 ,a), (15) 

Ri,R 2 >0,ae[0,l] 

where we have made explicit the dependence on (Ri, R 2 , a). As a benchmark, it is useful to 
consider the single-layer scheme obtained as a special case of the proposed strategy with R 2 = 
and a = 1. Thus, the optimal throughput of a single-layer strategy is the solution of the following 
problem: 

maximize f](Ri, 0, 1). (16) 
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IV. Constant Compression Gain 

In this section, we analyze the throughput of the proposed multi-layer HARQ strategy when 
the relay uses a constant compression gain a t = a Dt as in © for all t — 1, . . . , T regardless of 
the feedback information reported from the decoder. As explained in Sec. HU with this choice, 
the description W t can be recovered at the decoder for all realizations of the fading channel 
Sf However, this approach is not able to opportunistically leverage a more advantageous fading 
state St. A strategy that can exploit better fading state via adaptive compression will be discussed 
in Sec. |Vl We focus on the LTSC model, so that D t = D and S t = S for all t = 1, . . . , T. 
Moreover, we study both the case with local CSI at the encoder, i.e., when the encoder knows 
the local fading state D = d and thus can choose the tuple (Ri(d), R 2 (d), a(d)) as a function 
of d, and the case with no local CSI at the encoder. 

To express the objective throughput r] in (flOl) . we have to compute the probabilities in (fTJb 
and CHI) as a function of parameters R\, R 2 and a which is done in the following lemmas. 



6(d) 



Lemma 1. The probability pl ut (k) with compression gain cld is given as 

pl ut (k)=E D [6(D)} (17) 

where the function 9(d) is defined as 

M$§&+e)> if ci(^, i) < o 

1, if Ci(d,l) >o' 

with bd = 1 + dd/d and the functions Q(d, I) and Ci(d, I) given as 

C<(d,0=2 2fli(d)/, -l, (18) 
and (i(d, I) =2 2R ^ d)/l a(d) - 1, (19) 

for 2 = 1,2 and k — 1 , . . . , T. 

Proof: The proof is in Appendix \A-A\ 

m 

Lemma 2. If aP 1, the probability pl ut (k) with compression gain ao is approximated as 



(20) 



i=i 
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where v\ ,B {d) and v^J(d) are defined as 



oo, 



and <f(d)=minV fe u * 



ClW) _ OA 

Ci(d,l)P b d 



QB 



v fc u f = R u FW}u 

We have defined the function v^f(d) as 



CiW-i) _ ^ 



v V k f(d) 



oo, 
0. 



(C2(rf,fc-/)+l)^ /(fc -°- 1 a d 1 



Pc(d)'/(fc-0 b d P 

wjY/z the function c(d) given as c(d) = bd + a(d)Pa<i. 
Proof: See Appendix \A-B\ 



if Ci(d,Z) < o 
if Ci(<M) > 



if Ci(d,l- 1) > 
, ifCi(d,Z-l)<0 



it k = / and — d > 1 

if Jfe = Z and 22% !!^ < 1 



c(d) ; 



if k > I 



(21) 
(22) 



(23) 



(24) 



Lemma 3. If aP <C 1, ?/ze probability p\ cc {k) with compression gain an is approximated as 



k-l 



Pdec(fc) ~Ys Ed WW) - ^(«S(I>))]" 



+ E 



[^(g, UB P))-^(g fc LB P))] 



i=i 



where u^ B (d), u^f(d), q\ (d) and q^ B (d) are defined as 

u]^ B (d) =m&xUl 



/LB 

*k,l ' 



ttSJ(d)=minWjf, 

LB / a /-.LB 



Cx(MziL_M i fC 1 (rf,A;-l)<0 



and g fc UB (^) = ^ fl(d ' fc - 1)P 6d 

oo, if Ci(d, fc — 1) > 



(25) 
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U7 



ith the sets U^f, U^f and Q^ B given as 



{oo}, ifGKO>o 



CiW-i) 



. , „ 2d L if A(cU - 1) < 
0, if Ci(d,Z — 1) > 



and Q 



LB 



(2^, fc) _ Orf 

a(d)P bd 



U 



Ci(d,fe)P 6 d J ' 



{(X)}, 



Wfe /zave defined the functions u]^{d) and u^f (d) as 
(( 2 (^-l) + l)6f- lH 



Pc(d)J/(*-0 



6 d 



1 
P 



oo, 




2d _ _L 



Proof: The proof appears in Appendix IA-CI 



if &(d,k) < 
if Ci(d, k) > 



c(d) ! 



if / < jfe- 1 



(26) 



(27) 



(28) 



if Z = - 1 and — c ^ d y d > 1 
ifl = k — l and ^ < 1 . (29) 



With Lemmas [UGl we can express the throughput (flOt as a function of the optimization 
variables R2 and a via numerical integration over the distribution Fo{d). The optimization 
problems (TT3T ) and ([TBI are not convex and need to be solved via global optimization tool such 
as the branch-and-bound method Il20l . Specifically, with local CSI at the encoder, one needs to 
optimize over the parameters Ri(d), R,2{d) and a(d), which corresponds to the tuple (Ri, R2, a) 
to be used when the relay fading state is D t = d. In practice, this optimization can be reformulated 
by quantizing the fading distribution. Instead, without local CSI at the encoder, the optimization 
is done over the single tuple i?2, a) since the encoder is not able to adapt to the fading 
state D. 

V. Adaptive Compression Gain 

In the previous section, we have assumed that the relay employs Gaussian test channel © 
with compression gain a t = an for all t = 1,...,T regardless of the feedback information 
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reported from the decoder under the LTSC (|3). We recall that this choice guarantees reliable 
decompression even in the worst-case fading state S t = s min . This section is motivated by the 
attempt to leverage better fading states when they occur. To this end, we assume that the feedback 
information that only the message M% of the first layer was decoded in a slot k implies that the 
fading coefficient S t of the side information Y t is larger than some level s, that is, S > s. This 
can be calculated as 

f [{2 2R ^ k - l)b D + {2 2R ^ k a - l)Pa D \ 
| b D {\ -2 2R ^/ k a)P 

by imposing the condition that the accumulated mutual information is sufficient to support rate 
Ri (see Appendix |A] for more discussion). Upon reception of a positive acknowledgement for 
layer 1 and a negative acknowledgement for layer 2, we then propose that, from the next slot 
t > k, the relay performs compression assuming the better side information S t = s > s m i n . The 
corresponding compression gain is given as 



max 



(30) 



a D . 



(31) 



l/D + (l + s/D)P 

With adaptive compression, the expected throughput rj in (flOl) can be computed using the 
lemmas presented in the previous section with the changes discussed in the following lemmas. 

Lemma 4. If aP 1, the probability pl ut (k) with adaptive compression is approximated as 
([20]) in Lemma \2\ with v^J{d) modified as 



oo, 
0. 



^ 2 {d,k~l)+l)bf k - l) _ ~ ad 



Pc(d) i /(fc-0 



it k = I and — c ^ d y d > 1 

22R 2 (d) b l 

i£k = l and c{d)i d < 1 , 
if k > I 



(32) 



where bd = 1 + aa/d. 

Proof: See Appendix \B-A\ 



Lemma 5. If aP 1, the probability p^^k) with adaptive compression is approximated as 
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f |25l ) in Lemma\3\with u]^j(d) and u^f(d) modified as 



(( 2 (d,h-l) + l)b 



i/(k-i) 

d 



p c (d)i/(k-D 



1 

P 



(33) 



u V k f(d) 



oo, 



22R 2 (d) b l 

ii I = k -1 and — 7^-^ > 1 

c(dy 

if I = k - 1 and 2 £3 < 1 



c(d) ; 



(34) 



(C2(rf,fc-l-Q+l)6rf /(fc ~ 1 ~° 
Pc(d) ; /(fc-i-0 



if I < k - 1 



Proof: See Appendix \B-B\ 



With the results in Lemmas [Q 0] and |51 we can express the average throughput with the 
adaptive compression described in this section as a function of the design parameters R\, R 2 
and a. 

VI. Short-Term Static Channels 

In this section, we discuss the STSC model in which the channel coefficients D t and St change 
independently from block to block. For simplicity, as in IfPTTl . we focus on the case T = 2, i.e., 
there can be at most one retransmission. It is observed that, even with T = 2, we have to 
consider four random variables D\, D 2 , Si and S 2 , which complicate the analysis as compared 
to the LTSC model. Moreover, given the independence of the channel fading gains from block 
to block, adaptive compression is not applicable under the STSC model. Therefore, we set the 
compression gains as a t = ao t in © for t = 1, 2. The quantities in (ITT]) and CLU) reduce to 



E [R] = i2i (1 - pL t (2)) + R 2 (1 - pL(2)) 
andE[L]=pL(l) + 2(l-pL(l))- 



(35) 
(36) 



Thus, it is enough to compute three probabilities p l out {2), pi ut {2) and p 2 dec (l), which are derived 
in the following lemmas. 

Lemma 6. The probability p out (2) in the STSC model with T = 2 is given as 



pL t (2)=E DuD2>Sl [<?(£>!, £> 2 A)] 



(37) 
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where we have defined the function g(di, d 2 , Si) as 

) -l)b d2 +{2 2h ^ d i' B ^a 
(2 2h ( d l> s l)a-l)b d2 P 



Fq , [(2^^)-l)^ 2+ (^^-i)^-l)a d2 P]- . 22h(dl , si) j 



g{di,d 2 , si) = < 



if a(di, Si) = and a = 1 
otherwise 



with the function a(di,si) given as 

h(d 1 ,s 1 ) =Rx - fi (aP,aP,a dl ,si,di) . 

The function fi(P, P, a, s, d) is defined as 

f (p r> j\ 1 i (a , v> s + a(l + s/d) 

fi{P, P, a, s, d = - log 1 + P — 

2 \ 1 + a/a + P (s + a(l + s/d)) 

Proof: See Appendix \C\ 



(38) 



(39) 



(40) 



Lemma 7. The probability pl ut (2) in the STSC model with T = 2 is given as 

p 2 out (2) =E Dl , D2 , 5l [<p(D 1 ,D 2 ,S 1 )]+E DuD2! s 1 [ 1 (D 1 ,D 2 ,S 1 )]+p 1 out (2), 
where the functions tp(di, d 2 , s{) and 7(di, d 2 , s{) are defined as 

" 2 2(«2-/l(ai',0,a dl , S i.dl))_ 1 



ip(d 1 ,d 2 ,s 1 ) 




(41) 



and 7(di, d2, si) 



, if fi (aP, aP, a dl , s u d x ) > R 1 

if fi (aP, aP, a dl ,s 1 ,d 1 ) < R 1 

(42) 

[Fs 2 (a*ub(si)) - ^ Gulb(si))] + , if /i (a-P, aP, a dl , s u di) < i?i 



0. 



if /i (aP, aP, a dl , s ± , di) > i?i 



(43) 
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with /iub(<Si) and /xlb(si) given as 



/^ub(si) = < 



oo, 
0. 



if a = 0and2 2 ( fl2 - /l (" p '°^i' Sl ' d 0) > 1 
if 5 = 0and2 2 ( i?2 - /l (" p '°' ad i' Sl ' dl )) < 1 



if a > 



! 2(fl 1 -/ I (aP,aP,a dl , S1 ,d 1 ))_ 1 



b d2 -aPa d2 



1-2 v 1 : v d i 1 v/a )b d2 P 



and /xlb('S'i) = < 



0. 



\ 



(44) 



j£ 2 2 (^i _ /l( QP ' QP ' a di' s i' d l))Q, < I 



if 2 2(fi 1 -/ I ( a P,«P,a il , Sl ,d 1 ))- = 1 . 



and < aPddo 



otherwise 



Proof: See Appendix \C\ 

Lemma 8. The probability p\ ec iX) in the STSC model with T = 2 is given as 

p 2 dcc (l)=E Dl [max{A 1 p 1 ),A 2 (L> 1 )}] ) 



with the functions Ai(di) and A 2 ((ii) given as 
Aipi) = 



[(2^i-l) bdl+ (2^i«-l)Pa dl ] + 2B 

(l-2 2i? ia)b dl P ' 11 Z « <- J- 



OO, 



if 2 2Rl a > 1 



fed, 



and A 2 (Di) = < 



oo, 



0. 



if a = and i? 2 > • 
if a = and i? 2 = 



Proof: See Appendix \ 



(45) 



(46) 
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SNR p D of local fading D [dB] 

Figure 2. Average throughput -q versus the SNR p D with T = 2, C ma x = 1, P = OdB, p s = OdB and K = 0. 

VII. Numerical Results 

In this section, we present numerical results to gain insights into the advantage of the proposed 
multi-layer HARQ strategies. In the figures, the cases with and without local CSI at the encoder 
are denoted by "LCSIT" and "No LCSIT", respectively. We assume Rayleigh fading for the side 
information Y t and Rician fading for the signal V t received by the relay with Rician factor K (i.e., 
K is the ratio of the power of line-of- sight (LOS) component to that of non-LOS component). 
The rationale behind these distributions comes from the application to the cloud radio access 
scenario (see Sec. H]), in which Y t is the signal received by out-of-cell BSs, which typically 
lack the direct LOS component, unlike the signal V t received by the in-cell BS. The signal- 
to-noise ratios (SNRs) of the source-to-relay and the source-to-destination links are defined as 
E[|A| 2 ] = Pd andE[\S t \ 2 ] = p s , respectively. 

We first examine in Fig. |2]how the SNR pp of the relay fading channel D affects the average 
throughput rj by plotting rj versus the SNR p D under the LTSC model with T = 2, C max = 1, 
P = OdB, ps = OdB and K = 0. With local CSI at the encoder, the proposed BC scheme 
shows performance gain over the conventional single-layer approach only in the range of low 
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Figure 3. Average throughput r\ versus the backhaul capacity C max with T = 2, P = OdB, pr> — ps = OdB and K = 0. 

SNR p£>. This is because in this case, BC is only used to combat the uncertainty of the fading 
gain St, whose relevance becomes less pronounced as po increases. However, with no local CSI 
at the encoder, the gain of the BC remains substantial for all SNRs p D , since in this case, the 
CSI uncertainty at the encoder includes both D t and St. 

In Fig. |3l we plot the throughput performance versus the backhaul capacity C max for the LTSC 
model with T = 2, P = OdB, p D = p s = OdB and K = 0. It is observed that the impact of 
the local CSI at the encoder becomes more significant for larger backhaul capacities C max , since 
the performance is more affected by the encoder-to-relay link if the backhaul capacity C max is 
large enough. Moreover, the flexibility afforded by BC makes the effect of LCSIT less relevant 
than for conventional single-layer transmission. 

In Fig. HI we observe the effect of the maximum number T of transmissions for the LTSC 
model with C max = 1, P = OdB, p D = 10 dB, p s = OdB and K = 0. For both cases with 
local CSI at the encoder or not, the advantage of the BC scheme diminishes as T increases. 
This implies that the HARQ strategy is able to compensate for a large fraction of performance 
degradation of the single-layer scheme when enough number of transmissions are allowed. This 
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maximum number T of transmissions 

Figure 4. Average throughput r\ versus the maximum number T of transmissions with C max = 1, P = OdB, po = 10 dB, 
ps = OdB and K = 0. 

trend is more apparent in the case with no local CSI at the encoder, due to the layer gains of 
BC. 

In Fig. [51 we investigate the advantage of the adaptive compression scheme proposed in Sec. [V] 
by plotting the throughput performance versus the Rician factor K for the LTSC with C max = 2, 
T = 2, P = dB and p D = p s = 20 dB. We recall that the adaptive compression was proposed 
to opportunistically leverage better fading states. In accordance with this motivation, the adaptive 
compression is observed to be advantageous as the factor K grows due to increased frequency 
of good fading states that can be exploited via the proposed strategy. 

Finally, in Fig. [6l we compare the average throughput performance under the LTSC and STSC 
models with no local CSI at the encoder and T = 2, P = dB, C max = 5 and K = 0. With 
single-layer transmission, it is seen that the STSC model leads to better performance than LTSC 
due to the diversity gain. However, with BC, an additional factor determines the performance 
comparison, namely the possibility for "opportunistic retransmission" under LTSC. Specifically, 
under the LTSC model, when the encoder is reported an ACK for the first-layer message Mi, it 
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Figure 5. Average throughput r\ versus the factor K of Rician fading with C ma x = 2, T = 2, P = dB and po = ps = 20 dB. 




Figure 6. Average throughput r\ versus the SNRs po = ps with no local CSI at the encoder and T = 2, P = dB, C" m ax = 5 
and K = 0. 
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can transmit the second layer M 2 in the next slot in order to leverage the good fading state. In 
contrast, under the STSC model, this is not possible since the fading coefficients D t and S t vary 
independently from block to block. From the figure, it is observed that this factor is dominant 
in the low-to-moderate SNR range, where the performance of BC transmission under the LTSC 
model is better than under the STSC model. 

VIII. Conclusions 

Motivated by the uplink of cloud radio access networks, we have studied robust transmission 
and compression schemes for the fading relay channel with an out-of-band relay. Specifically, 
we have adopted a multi-layer BC transmission strategy coupled with HARQ, thus allowing 
for variable-length transmission and variable-rate decoding, under two different channel models, 
LTSC and STSC. Moreover, we have proposed an adaptive compression strategy at the relay that 
is able to leverage better fading state based on the HARQ feedback received from the destination. 
We have demonstrated the performance gain of the proposed schemes over conventional single- 
layer approaches via extensive simulations. 

Appendix A 
Proof of Lemmas |T}{3] 

In this Appendix, we derive the probabilities presented in Lemmas [Q13] Since we assumed the 
LTSC model © in Sec. El we have D t = D and S t = S for alH = 1, . . . , T. We first calculate 
the probabilities conditioned on D = d and the results in Lemmas [D|3] are then obtained by 
taking expectation over D. 

A. Proof of Lemma [7] 

In this subsection, we compute the probability pl ut (k) that the message M\ is not decoded 
until slot k. Since we have assumed the IR-based HARQ approach, the probability pl ut (k) can 
be calculated as the probability that the mutual information accumulated along the first k slots 
is smaller than i?i lfT3l : 

pl nt (k) = E D [Pr [k ■ I(X ltt ; F t |BCmode) < Ri(D)\D]] , (47) 
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Pout 



1=1 



Pr 



where the mutual information I(X l t ;Y t \BCmode) is given as 

/(Xi, t ;y t |BCmode) = fi(a(D)P,a(D)P,a D , S, D), (48) 

with the function fi(P, P,a, s,d) defined in (1401) . If we express the probability (|4VT > using the 
CDF of S, we arrive at the expression (fTTt . 

5. Proof of Lemma [2] 

This subsection computes the probability p 2 OVLt (k) that the message M 2 is not decoded until 
slot k. Using the total probability theorem, we can write the probability p 2 ut (k) as 

(I - 1) ■ I(X ht ;Y t \BCmode) < R^D) < I ■ I(X ht ; Y t \BC mode), 
^ Z ■ I(X 2jt ; F t |X ljt , BCmode) + (k — I) ■ I(X 2it ; Y t \X u , SLmode) < R 2 

+ Plt(k), (49) 

where the probability pl ut (k) was derived in the previous subsection and the mutual information 
quantities related to the second layer are given as 

I(X 2>t ;Y t \X 1;t , BCmode) =f 1 (a(D)P,aP,a D , S, D), (50) 

I(X 2>t ;Y t \X 1>t , SLmode) =f 1 (P, 0, a D , S, D). (51) 

The term inside the summation in (|49l is then derived as 

(/ - 1) ■ Ipf ljt ;F t |BCmode) < Ri(D) < I • J(X ljt ; Y t \ BCmode), 
I ■ I(X 2 , t ; Y t \X 1>t , BC mode) + (k — I) ■ I(X 2it ; Y t \X 1>t , SL mode) < R 2 {D) 

(1 _ 2 2R ^ D y^a(D)) b D PS < (2 2 - R i( D )/C- 1 ) - l)b D + (2 2R ^ D y^a(D) - l)Pa D 
(2 2ij l( D)/i _ + (2 2 «iW/'a(D) - l)Pa D < (1 - 2 2R ^)/'a( J D))6 D P 1 S, 



Pr 



D 



Pr 



(1 + a(D)PS) + a(D)Pa D y (b D (l + PS) + Pa D f- 1 < 2 2R ^b k D 



D 



D 



(52) 



where the last condition makes it difficult to express the probability in terms of the CDF of S. 
Following IfTTll . we assume the low SNR condition a(D)P <C 1 so that we have l + a(D)PS ~ 1 
in the last condition of the probability (l52l . Then, the probability (l52l) is approximated as 

(Z - 1) -/(X^F^ BCmode) < R\(D) < I ■ I(X ht ;Y t \ BCmode), 
/ ■ 7(X 2)t ; Y t \X lit , BCmode) + (k — I) ■ /(X 2 , t ; Y t \X 1>t , SL mode) < R 2 (D) 

[Fs(v^f(D))-Fs(v^(D))] + , (53) 



Pr 



D 
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where v^f(D) and vffi(D) are defined in Lemma |2] If we substitute (1531 ) into (1491 ), the result 
in Lemma [2] is obtained. 



C. Proof of Lemma \3\ 

In this subsection, we compute the probability p\ ec (k) that the message M 2 is successfully 
decoded in slot k. Following similar arguments as above, we can write Pdec(^) as 



PdecO) 



k-1 



D 



1=1 



Pr 



(I - 1) ■ I(X u ;Y t \BCmode) <R X {D)<1- Ip^;^ BCmode), 
I ■ I(X 2tt ; Y t \X u , BCmode) + (k — 1 — I) ■ I(X 2/ , Y t \X 1:U SL mode) < R 2 (D), 
R 2 (D) < I ■ I{X 2>t ; Y t \X ljt , BCmode) + (k — I) ■ I(X 2>t ; Y t \X ljt , SLmode) 

(54) 



+ E 



D 



Pr 



{k-1) ■ /(Xi^^tlBCmode) < R^D) < k ■ I{X 1;t ;Y t \BCmode), 
R 2 (D) < k ■ I(X ljf ; F t |BC mode) 



D 



(55) 



Moreover, under the low SNR condition aP <C 1, the term inside the summation is approximated 



as 



Pr 



D 



(56) 



(I - 1) •/(Xi ) t;f t |BCmode) < R X {D) < I ■ I{X 1)t ;Y t \ BCmode), 
I ■ I(X 2tt ; Y t \X hu BC mode) + (k — 1 — I) ■ J(X 2 , t ; Y t \X u , SL mode) < R 2 (D), 
R 2 (D) < I ■ I(X 2it ; Y t \X hu BCmode) + (k — I) • I(X 2}t ; Y t \X hu SLmode) 

^[F s (u^(D))-F s (ul B (D))} + , 

where we have defined u)^f(D) and u^ B (D) in Lemma [3] Moreover, we can derive the last term 
in ([55]) as 

r (A;-l)-/(X M ;Fi|BCmode) < R X {D) < k ■ I{X ljt ] Y t \BC mode), 
R.2(D) < k ■ I(X ljt ; F t |BCmode) 

= [F s ( q r(D))-F s (q^(D))]\ (57) 

with q^ B (D) and q\ B (D) defined in Lemma |3j As a result, we obtain (|25l ) by plugging (1561 ) and 
(1571) into ([55]). 



Pr 



D 
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Appendix B 
Proof of Lemmas [4] and [5] 

In this appendix, we derive the results in Lemmas HI- |5] with adaptive compression. 
A. Proof of Lemma |?] 

If we assume the adaptive compression described in Sec. [V] the probability Po Ut (k) is calculated 
as <|49l with the mutual information I{X 2jt \ Y t \Xi >t , SLmode) changed from (|5ll to 

I(X 2>t ; Y t \X 1>t , SL mode) = A(P, 0, a D , S, D). (58) 

The only difference from (1511 is the improved compression gain old- If we calculate (|49l ) with 
(1581) . we immediately obtain the result in Lemma |4] 



B. Proof of Lemma \5\ 

With the adaptive compression, the probability p 2 dcc (k) is given as (1551 ), and similar to the 
previous subsection, the only difference is that the mutual information I(X 2jt ; Yt\Xi t , SLmode) 
is computed as (l58l) with the compression gain old- Then, we can obtain the result in Lemma |5] 

Appendix C 
Proof of Lemmas [MS] 

In this appendix, we avoid repetition by focusing on the proof of ([37]) in Lemma [6] since the 
proof for Lemmas [7][8] follows similarly. With STSC, the probability Pout (2) is given as 

Pout (2) 

= Pr [l(Xix, Yi|BG mode) + J(X 1>2 ; F 2 |BC mode) < i2j 

f 1 (a(D)P,a(D)P,a Dl ,S l ,D 1 )+ 
f l {a{D)P,a{D)P,a Da ,S 2 ,D 2 )<R 1 (D) 
If we express the conditional probability inside the expectation in (l59l) with respect to the CDF 
of the fading coefficient S 2 , we get Eq (1371) in Lemma [6l 



-El?i,Da,Si 



Pr 



D u D 2 ,S l 



(59) 
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